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Abstract 

Following an approach developed by Gieseker, Knorrer and Trubowitz 
for discretized Schrodinger operators, we study the spectral theory of 
Harper operators in dimension two and one, as a discretized model 
of magnetic Laplacians, from the point of view of algebraic geometry. 
We describe the geometry of an associated family of Bloch varieties 
and compute their density of states. Finally, we also compute some 
spectral functions based on the density of states. 

We discuss the difference between the cases with rational or ir- 
rational parameters: for the two dimensional Harper operator, the 
compactification of the Bloch variety is an ordinary variety in the ra- 
tional case and an ind-pro-variety in the irrational case. This gives 
rise, at the algebro-geometric level of Bloch varieties, to a phenomenon 
similar to the Hofstadter butterfly in the spectral theory. In dimen- 
sion two, the density of states can be expressed in terms of period 
integrals over Fermi curves, where the resulting elliptic integrals are 
independent of the parameters. 

In dimension one, for the almost Mathieu operator, with a similar 
argument we find the usual dependence of the spectral density on the 
parameter, which gives rise to the well known Hofstadter butterfly 
picture. 
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1 Introduction 



In mathematical physics, the Hamiltonian of a lattice electron in a uniform 
magnetic field, called the discrete magnetic Laplacian, has been studied for 
years in the tight-binding model [H]. As a special discrete magnetic Lapla- 
cian, the Harper operator corresponds to a square lattice when the coupling 
constant is fixed (i.e. A = 1) [6]. For instance, the Harper operator arises 
in the study of the integer quantum Hall effect [3]. One important spec- 
tral property of the Harper operator is that its spectrum is a Cantor set of 
zero Lebesgue measure for every irrational frequency [1], |12]. In this paper 
we show that the different structure of the spectrum (band or Cantor-like) 
in the rational and irrational case has an analog in terms of the algebro- 
geometric properties of the Bloch variety, which is an ordinary variety in the 
rational case and an ind-pro-variety (with properties analogous to a totally 
disconnected space) in the irrational case. 

More precisely, we consider the algebro-geometric approach to the spec- 
tral theory of electrons in solids developed for discrete periodic Schrodinger 
operators by Gieseker, Knorrer and Trubowitz in [5] and we investigate how 
to extend this approach to the case of Harper operators and almost Mathieu 
operators. As we will see in detail in ^below, the main obstacle, which so far 
prevented people from extending this method to the case with magnetic field, 
lies in the fact that, in the irrational case, one leaves the world of algebraic 
varieties, due to the simultaneous presence of infinitely many components. 
We will resolve this problem here by describing the relevant geometric space 
via a suitable "limit procedure" , by which one can make sense of spaces that 
are not themselves algebraic varieties, but that are built up from a sequence 
of algebraic varieties. Making sense of this limit procedure requires working 
in a setting where operations such as taking limits of geometric objects are 
possible, and this is the reason why we have to resort to a moderate use of 
categorical tools and methods, which are precisely built in order to allow for 
notions of limits of the type we need. 

Our goal in this paper is to obtain explicit expressions for the spectral 
densities of the two dimensional Harper operator and of its "degeneration" , 
the one- dimensional almost Mathieu operator, in terms of periods on a suit- 
able "limit of algebraic varieties" . In order to obtain this result, we first need 
to overcome two main problems of a geometric nature: 

1. Take a limit of algebraic varieties; 
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2. Compactify and remove singularities in a way that is compatible with 
this limit procedure. 

The first problem is the one we already mentioned: we will see that, 
in the case of irrational parameters, the geometric space that describes the 
spectral problem for the Harper operator is the locus of zeros of a countable 
family of polynomial equations. As such it is not an algebraic variety, but 
we will show that it can be constructed by taking a limit (in a suitable sense, 
technically called an ind-variety) of a family of algebraic varieties. Roughly 
speaking, this means that the relevant space can be built by progressively 
adding more and more components, so that one can carry out operations such 
as period computations essentially "component-wise", with certain compat- 
ibilities taken into account. 

The second problem is more subtle. In the original formulation of [5J , for 
the case without magnetic field, one needs to compactify the Bloch variety 
(which describes the spectral problem for the discretized periodic Schrodinger 
operator) by adding a Fermi curve "at infinity" and to resolve singularities 
by blowups, in order to be able to work in a smooth setting, where the pe- 
riod integral giving the spectral density function is computed. This problem 
occurs in exactly the same way when we extend this method to Harper and 
almost Mathieu operators, but it is now further complicated by the fact that 
we want to be able to perform these operations of compactification and re- 
moval of singularities in a way that is compatible with the "limit procedure" 
already mentioned, so that we can "pass to the limit" and still achieve the 
desired result. 

We will see that this problem forces us to replace a limit with a double 
limit. More precisely, on the one hand there is what one refers to as an 
inductive limit, namely the operation of assembling together the countably 
many components of the Bloch ind-variety. On the other hand there is the 
procedure that removes singularities in each component (again, with com- 
patibilities) and that is achieved by blowup operations. These create varieties 
that are progressively less singular and which project down to the original 
singular variety, in such a way as to not change anything where the latter was 
already smooth. This other operation defines what is called a projective limit 
and the desired compatibility of these two limiting operations that allow us to 
solve both problems (1) and (2) at the same time is encoded in what is called 
an ind-pro-variety in the current algebro-geometric terminology. Although 
we work in a very special case of this very general kind of construction, it 
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will be convenient to adopt this general framework, as that assures us that 
we can "pass to the limit" by consistently carrying out calculations in the 
intermediate steps, at the level of the successive approximations. 

While this geometric construction is appropriately formulated in algebro- 
geometric and categorical language where we can appropriately define the 
necessary double limit operation, the resulting space we obtain has a simpler 
heuristic description in more directly physical terms. In the original case of 
[5] of the periodic Schrodinger operator without magnetic field, the Bloch va- 
riety physically describes the complex energy-crystal momentum dispersion 
relation, that is, the set of complex points that can be reached via analytic 
continuation from the band functions. In the case of Harper and almost 
Mathieu operators with irrational parameters, one observes the phenomenon 
that, instead of intervals (band structure), the spectrum takes the form of 
a Cantor set, which is geometrically a suitable limit of an approximating 
family intervals. Correspondingly, the geometric space describing the com- 
plex energy-crystal momentum dispersion relation, which replaces the Bloch 
variety, is no longer an algebraic variety, but it is obtained in the same way 
by "analytic continuation" from a Cantor set rather than from bands in the 
spectrum. It is then natural to expect that such a geometric space will be 
the "algebro-geometric analog of a Cantor set", also obtained as a limit, in 
a way that mirrors the construction of the Cantor set as a limit of intervals. 
It is then amusing to see that this concrete finding is in complete agreement 
with a general categorical result of Kapranov [15], which shows that all 
ind-pro- varieties (which simply means all spaces that are obtained by a "dou- 
ble limit procedure" of the kind described above) in fact really do look like 
Cantor sets. 

The paper is structured as follows: the geometric problems (1) and (2) 
above are dealt with in §2.1| -^ |T51 A more physics oriented reader who is 
willing to believe that a suitable double limit of varieties needed to solve both 
the first and the second problem can be constructed, can skip directly to §2.61 
and the following sections, where the space obtained in the previous part of 
the paper is used to perform the needed period computations that describe 
the density of states and the spectral functions for the Harper and almost 
Mathieu operators. In particular, §2.61 and §2.71 deal with the case of the 
two-dimensional Harper operator, while ^ deals with the one- dimensional 
case of the almost Mathieu operator. 
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1.1 The algebro- geometric setup 



In [5] , Gieseker, Knorrer and Trubowitz modeled the theory of electron prop- 
agation in solids by considering a discrete periodic Schrodinger operator 
—A + V acting on the Hilbert space £^(Z^), with 

Atplm, n) = tplm + 1, n) + iplm — l,n) + ?/'(m, n + 1) + iplm, n — 1) (1) 

the random walk operator that discretizes the Laplacian, and with an effec- 
tive potential V given by a real function, periodic with respect to a sublattice 
aZ + 6Z, where a and b are distinct odd primes. They studied the geometry 
of the associated Bloch variety 



One then considers the projection tt : B{V) — t- C on the third coordinate 
and one defines the Fermi curves as Fx{C) := 7r"^(A). 

In more physical terms, the Bloch variety represents the complex energy- 
crystal momentum dispersion relation, that is, the locus of all complex points 
that can be reached by analytic continuation from any band function. 

Defining the integrated density of states p(A) as the averaged counting 
function of the eigenvalues, they observed that the density of states dp/d\ 
can be expressed as a period integral over the homology class Fx, namely 



where ux is a holomorphic form on the complex curve Fx{C). 

The interested reader is referred to 3D VRML Fermi Surface Database 
[1] for concrete physical examples of Fermi curves. 

1.2 The case with magnetic field: Harper operator 

In this paper we consider a different, but closely related, spectral problem, 
where instead of the usual Laplacian and its discretization given by the ran- 
dom walk operator, one considers a magnetic Laplacian, whose discretization 
is a Harper operator. We also restrict our attention to the case with trivial 
potential V = 0. 



B{V) 



{(a, 6, A) G C* X C* X C| 

3 nontrivial ip with (—A + V)ip = Xip, 

such that ip^m + a,n) = ^iip{m, n) 

and '?/'(m, n + h) = C,2'ip{^, n), V (m, n) G Z^ 



(2) 



} 
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The (two dimensional) Harper operator H acting on is defined as 



g-27ri/3m^(^^^ Tl + 1) + e2'^*'^'^'i/>(m, U - 1) 



where the two unitaries 



f/^(m, n) := e-2™"V'(^ + 1, n) and ■1/V^(m, n) := e-^'''^"'ip{m, n + 1) 

are the so-called magnetic translation operators with phases a and (3 re- 
spectively, and the group of magnetic translations T^^js is generated by U, V. 
One can then write the Harper operator as H = U + U* + V + V* . Note 
that our form of the Harper operator is slightly different from that in the 
literature, but they are all unitary equivalent by a gauge transformation 
T^ilj{m,n) = e2'^'T"""V'(m,n) on i'^{Z^). 

Let Ta,i3 = C*{Ta^i3) be the group C*-algebra of the group of magnetic 
translations Tq,^^. Recall that the noncommutative torus Ag is the universal 
C*-algebra generated by two unitaries u, v subject to the commutation rela- 
tion uv = e^'^'^^vu. Setting 9 = a — (3, we have a representation tiq : Ag ^ Ta,/3 
such that 7r6i(-u) = U,TCg{v) = V. Thus, the Harper operator H is the image 
of a bounded self-adjoint element of Ag. 



1.3 Rational vs irrational 

If the parameters a, /3 are rational numbers, so is 6, then the rotation algebra 
Ag is isomorphic to the continuous sections of some vector bundle over the 
two-torus T^. For rational parameter, the one dimensional Harper operator 
is a periodic operator and its spectrum consists of energy bands separated 
by gaps by Bloch-Floquet theory. The Bloch variety associated to a Harper 
operator with rational phases has only finitely many components, and can 
be treated similarly to the periodic Schrodinger operators discussed in [5]. 

When 6 is an irrational real number, the irrational rotation algebra Ag 
is a simple C*-algebra and it has been studied in noncommutative geometry 
motivated by Kronecker foliation, deformation theory etc. For irrational pa- 
rameter, the spectrum of the Harper operator is a Cantor set of zero measure 
[6], 0. 

Here, instead of following the approach to the spectral theory of Harper 
and almost Mathieu operators based on functional analysis and noncommu- 
tative geometry, we aim at adapting the algebro-geometric setting developed 
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in [5] in the case of discretized periodic Schrodinger operators without mag- 
netic field. 

As we show in [|2] below, in the case of irrational parameters, the analog of 
the Bloch variety that describes the complex energy-crystal momentum dis- 
persion relation is no longer an ordinary algebraic variety, since it is defined 
by a countable family of polynomial equations. It can still be described in 
such a way that algebro-geometric methods apply, as a "limit" of algebraic 
varieties, an ind-variety that has infinitely many components. This descrip- 
tion brings in some additional difficulty when we compactify it and blow it 
up to resolve the singularities. 

This second operation requires one more limit operation, which needs to 
be compatible with the first one. An algebro-geometric setting where such 
double limit operations can be described and the necessary compatibilities 
are encoded in the structure is that of ind-pro-varieties, which can be viewed 
as locally compact Hausdorff totally disconnected spaces, and behave in all 
ways more like Cantor sets than algebraic varieties (see [9], [I5]). Thus, this 
different behavior reflects at the level of Bloch varieties the different structure 
(bands or Cantor sets) of the Hofstadter butterfly spectrum in the rational 
and irrational cases. 

We assume that a, /3, 9 are all irrational real numbers in this paper, which 
will not be stated otherwise. 

1.4 Harper and almost Mathieu operators 

A limit case of the Harper operator is the almost Mathieu operator, whose 
spectral theory has been widely studied in connection with the famous phe- 
nomenon of the Hofstadter butterfly, which was flrst observed by Hofstadter [6], 
also cf. 0, [13]. 

Indeed, if we set the parameter /3 = and let the Harper operator act on 
£^(Z), we can express the resulting operator in terms of two new unitaries 

U'^in) ■= e-2™Xn) and V'^{n) := ^{n + 1) 

Thus the almost Mathieu operator is deflned as iJ' := f/' + U'* + V' + V'*, 
namely 

H'ip{n) := 2cos{27Tan)ip{n) + ip{n + 1) + ^{n - 1) (4) 

The physical meaning of this limit process is that one takes the Landau gauge 
which forces the vector potential only in one direction, so the almost Mathieu 
operator is sometimes called the Landau Hamiltonian by physicists. 
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Therefore, we get another representation tTq, of the noncommutative torus 
such that Tiaiu) = U',TTa{v) = V . In the hterature, the almost Mathieu 
operator is also referred to as the (one dimensional) Harper operator. Since 
the C*-algebra Ag is simple for irrational 6, the almost Mathieu operator has 
the same spectrum as that of the (two dimensional) Harper operator. 

The spectral property of the Harper operator is related to the famous 
"Ten Martini problem", which was one of Simon's problems ^7\. Now its 
spectrum is well understood and the reader can consult the sources listed 
as follows. Applying semi-classical analysis based on Wilkinson's renormal- 
ization, Hellfer and Sjostrand [8] proved that the spectrum of the almost 
Mathieu operator has a Cantor structure for frequency a having continued 
fraction expansions with big denominators. Last [T?| also showed that it is 
a zero measure Cantor set for a satisfying a Diophantine condition. 



2r 1 1 r 




p/q= 

Figure 1: Energy spectrum (—2 < < 2) as a function of the magnetic flux 
a = $/$o where rational numbers < p/q < 1 are used to approximate a. 

1.5 Spectral density and parameters 

In the discrete approximation, we study the geometry of a family of Bloch 
varieties and Fermi curves associated to the spectral theory of the two di- 
mensional and the one dimensional Harper operators, and we compute their 
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density of states. For the (two dimensional) Harper operator, we show in 
§2.61 that the density of states turns out to be a period independent of the 
parameter 9 = a — P, which involves the complete elliptic integral of the 
first kind. However, for the almost Mathieu operator (one dimensional case), 
we show by a similar technique in ^ that the density of states explicitly 
depends on the parameter a in a way that recovers the usual phenomenon 
of the Hofstadter butterfiy, illustrated in the figure. 

Note that, besides the difference in dimension, the generators of the 
Harper operator have a symmetric form in the variables and parameters, 
while the almost Mathieu operator is naturally generated by translation and 
rotation operators. 

Because of the electron-hole symmetry, the density of states of the Harper 
operator takes the same form for positive and negative energy levels. We 
show in §2.61 that it can also be represented as a sum of half-periods of two 
isomorphic elliptic curves. 

When the density of states is available, some spectral functions, such as 
the zeta function and the partition function, can be obtained by integrations 
of special functions of the eigenvalues. 

2 Harper operator and algebraic varieties 

In this section, we first describe the geometry of a family of Bloch varieties 
and algebraic Fermi curves associated to the Harper operator in the discrete 
model. Our discussion here follows closely the setting and notations used in 
[5] and we will point out explicitly where and how our case differs from the 
case without magnetic field considered there. 

While the rational case is essentially like the original case of [S], in the 
case of irrational parameters new phenomena arise, which refiect the different 
structure of the spectrum (Cantor sets instead of bands). 

We first show that the analog of the Bloch variety, which describes the 
complex energy-crystal momentum dispersion relation, is no longer defined 
by a finite set of polynomial equations but by a countable family of such 
equations. Thus, it is no longer an algebraic variety, but we can still describe 
it as an inductive limit of algebraic varieties (an ind-variety) by a procedure 
that essentially amounts to working with an approximating family of alge- 
braic varieties obtained by considering only finitely many components at a 
time, with compatibility conditions. 
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We then deal with the problem of the compactification and singularities. 
Using similar embedding and blowup techniques as in §5 and §7 of we de- 
scribe the compactification of our Bloch ind- variety. This requires performing 
the necessary operations compatibly on the successive approximations of the 
ind-variety. The resulting system of operations can be described as a double 
limit procedure and it is encoded in a compatible system of embeddings and 
projections. This is a special case of a general procedure devised precisely 
to the purpose of making sense of such geometric double limits, which is 
called an ind-pro- variety, and which in our case gives the algebro-geometric 
counterpart of the Cantor sets in the spectrum, in the form of one of these 
ind- pro-objects that, although obtained from families of algebraic varieties, 
behave themselves like Cantor sets. 

We carry out this construction in the subsections §2.1| - §2.5[ We then 
focus, in the later subsections starting with §2.6[ on the computation of the 
density of states and results based on that, such as the spectral functions, 
using period integrals computed on the spaces constructed in §2.1| - §2.5[ 

2.1 A family of Bloch varieties 

As in [5j , one defines the Bloch variety associated to the Harper operator ([3]) 

as 

B := {(6,6,A) G e X C* X C|i/^ = AV^, 

ip{m + a,n) = ^iijj{m,n), ^{m,n + b) = 

This looks very similar to the original case ([2]) with trivial potential V = 0. 
However, in fact, the locus determined by (EJ differs significantly from the 
case without magnetic field, since in the irrational case it defines a countable 
collection of algebraic varieties. 

To see this, we first describe the spectral problem defining S as a count- 
able collection of matrices acting on the vector {il){m,n)} e C"'', when a 
and /3 are irrational. This becomes a finite family when both a and /3 are 
rational. 

The correct notion of "limit of algebraic varieties" that we need to employ 
here to describe the geometric properties of the resulting space is that of an 
ind-variety. We first recall the abstract definition and then we give a more 
heuristic description of its meaning, before we apply it to our concrete and 
specific problem. 
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Definition 2.1. An ind-variety over C is a set X together with a filtration 

X° C C C . . . 

such that U„>oX" = X and each X" is a finite dimensional algebraic variety 
with the inclusion X" X""*"^ being a closed embedding. Such ind-variety 
X will also be denoted by limX". An ind-variety X = lii^X"' is said to be 
affine (resp. projective) if each X" is affine (resp. projective). 

In categorical terms, an ind-variety X is a formal filtered colimit of an 
inductive system {X"} of varieties. The notion of an ind-variety was first 
introduced by Shafarevich and one can find more examples and properties of 
ind- varieties for instance in |lUj . 

Let X and Y be two ind- varieties with filtrations {X^} and {i^"}. A 
map / : X — i- y is a morphism if and only if for every n > 0, there exists 
a number m(n) > such that /Ix" : X"^ — > is a morphism between 

varieties. Thus we get the Ind category Ind(V^ ar) of the category of varieties. 

Clearly, the idea here is that one "builds up" the limit space by a sequence 
of step-by-step operations that progressively enlarge an algebraic variety. 
The limit space is no longer a variety itself, but it can be approximated by 
varieties and one can carry over algebro-geometric operations to the limit, as 
long as they are performed on the approximating varieties in a way that is 
compatible with the successive inclusions that build up the limit space. For 
example (as will be directly relevant to our case) the limit may be obtained by 
adding more and more components. The abstract definition recalled above is 
much more general and allows for more complicated cases where, for example, 
the dimensions of the successive approximations grow, so that the limit may 
also be infinite dimensional. Our case is milder, as the dimensions of the 
components remain bounded and only their number is growing. However, it 
is convenient for us to describe the result as an ind-variety as that encodes all 
the compatibility conditions between the successive approximations that we 
will need to respect when dealing with the compactification and singularities 
problem. 

We now proceed to see what all this means concretely, in our specific case. 

Lemma 2.2. The Block variety ^ defined by the spectral problem of the 
Harper operator is an affine ind-variety that can be written as B = 
Uk,e(^zBk,e, where 

Bk,i = {(6, 6, A) G C* X C* X C I Pk,M,, 6, A) = 0}, (6) 
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for polynomials 



PM(ei,6,A) = det(M(^'^)-A/), 



(7) 



where the ah by ah matrix M*^'^'^^ has entries 



if m' = m,n' = n 




(8) 



Proof. Consider the spectral problem Hil){m,n) = \ip{m,n). Since we work 
with the boundary conditions 'il){m + a, n) = ^iip{m,n) and ip{m,n + h) = 
^2'ip{nT', n), we can consider the range where m = 1, . . . , a and n = 1, . . . ,h. 
If a and /3 are irrational numbers, then the phase factors exp{27iian) and 
exp(27ri/9m) are not periodic. This means that, for each (m, n) in the chosen 
fundamental domain of the aZ © 6Z action, we have a collection of problems, 
parameterized by the choice of an element {k,i) in Z^, which differ only in 
the presence of the phase factors exp(27ria(r;, + ih)) and exp(27ri/3(m + ka)). 
In the case where a is rational, the phase factors exp{27Tia{n + ih)) repeat 
periodically, with only finitely many distinct values, and so for exp{27i i/3{m + 
ka)), when /3 is rational. Thus, in the case where both a and (3 are rational, 
there are only finitely many different varieties B^ i to consider, so their union 
is a genuine algebraic variety. While in the case where at least one of the two 
parameters is irrational there are infinitely many components, so their union 
B is an ind-variety. For each pair {k, i) then one can write the corresponding 
problem in the form given by the matrix ([8]), by arguing as in §2 of [5]. In 
addition, by the form of the defining matrix, B^ i is symmetric under the 
involution on each fiber, A) (-)■ {^i^, A). 



The results of this subsection take care of the first of the two geomet- 
ric problems mentioned in the introduction. In the following subsections 
§2.2| - §2.5I we deal with the second problem: compactification and singular- 
ities. We show that we can resolve it in a way that is compatible with the 
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limit procedure described in this subsection. This will be the more technical 
algebro-geometric part of the paper, and it can be skipped by the readers 
who wish to see directly the derivation of the density of states and spectral 
functions in terms of periods, and who are willing to assume, without going 
through a more detailed explanation, that we can indeed make sense of the 
double limit procedure needed to obtain the correct geometric space on which 
these period computations take place. 



2.2 Fourier modes description 

In this subsection we just rewrite our varieties with a convenient change of 
coordinates, coming from Fourier transform, which will be useful later, when 
we deal with the compactification and singularities problem. 

Let the Fourier transform of tfj e be G L^(]R^/Z^), namely 



m n 



where Zi := e'^'"^^^,Z2 := e^'^*'^^ and a^n '■= ipim^n). Furthermore, the Fourier 
transform of Hip is 

+e-^^'^^^{ki - 13, k2) + e^^'^^^i^ki + /3, fca) ' 



I.e. 



HiP{zi, Z2) = Yl (e"^™"^r^ + e^'^^""^! + e-2"^^™22"^ + e'^—Z2)amnzrz2 



(10) 

Thus the spectrum of H is given by the Fourier modes in the fc-space: 
If we consider functions in £^(Z^/aZ + bZ), then by introducing {k,i) G 

^ _ ^-2iTia{n+eb) _|_ ^2iTia{n+eb) _|_ g-27rj/3(m+A:a) ^-1 _|_ ^2iTil3{m+ka) ^j^j^-j 

We then proceed, for a fixed {k,i), as in §2 of [5]. We introduce the 
unramified covering 

c : C* X C* X C ^ C* X C* X C 
{zi,Z2,X) ^ {z1,zlX) 
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and the preimage B^/ := c~^{Bkg). So the covering c : Bj^^g — )■ B^i has 
the structure group fia x fJ'b, where by /i„ we mean the group of roots of 
unity of order n, with the action of p E fia x l^b on the fibers of the form 
p- (zi, Z2, A) = {piZi, P2Z2, A). We write P = P o c. 

We now show how the discrete Fourier transform description of the Bloch 
varieties given in §2 of [5] is affected by the presence, in the Harper operator, 
of the phase factors exp(27ria) and exp(27r2/3) and their powers. 

For {m,n) G Z^, {z'^zl^} consists of a basis for the functions in L^(T^). 
There is an obvious action of pa x Pb on this basis p ■ z'^z'^ := p^p^z^z"^, 
which is basically a change of base, and can be naturally extended to an 
action on ^^(T^). 

Let p act on HiIj{zi, Z2). It not only changes the basis from {z'^zl^} to 
{pTP2^T^2}y but it also changes the Fourier modes into 

Fix {poi, P02) = (e^'^*/", e^'''/*), other roots of pa and p^ can be written as 
{pi, P2) = (Poi)Po2) some integers 1 < p < a and 1 < q <b, then rewrite 
the Fourier modes as 

-a{n+eb)a-p I a(n+£b)a+p -/3(m+fca)fe-g 1 8{m+ka)b+q 

Poi ^1 + Poi ^1 + P02 ^2 + P02 ^2 

Lemma 2.3. The Harper operator ([3]) determines a family of operators 
}j{k/) ^ /or G 7?, which acts as multiplication by the complex number 

a{n+eb)a+p -a{n+eb)a~p I I3(m+ka)b+q ~l3{m+ka)b-q I /^r)\ 

POI ^1 +P02 ^2+Po2 ^2 U^J 

Proof. In addition to the above discussion, we also have to take care of the 
boundary conditions. Let us look at one of them il){m + a, n) = C,i'ip{m,n). 
Taking the Fourier transform on both sides gives e~^'^*°'^^^/'(A;i, ^2) = C,iip{ki, ^2) 
or equivalently ip{zi,Z2) = z^^iip^zi, Z2). So the boundary conditions are re- 
moved naturally in the covering space Bk,e if we set ^1 = z^"" and ^2 = Z2^- 
In fact, take the involution symmetry (^1, ^2, A) 1— )■ (^f^, ^2^^, A) into account, 
we can just set ^1 = z^ and ^2 = ^2. □ 

We write M*^*^'^) for the diagonal ab x ab matrix with entries 

M^i\p, z) := p^J"+'')"+^;^i+PoT^"^''^"""^r'+Pof'''"^'^'^2+Pot^'"^''^^''"'^2"' 

(13) 
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Thus for pairs {k, i), B^/ is the zero-set of Pk/{zi, Z2, A) = det{M^^'^^ — A/), 
or equivalently it is given by the zero locus 

Bu,, = {{z,,z,,X) I n (M^iKp,^) - A) = 0}. (14) 

m,n,p 

For brevity, hereafter we denote g := ((m,n), {pi, P2)) G x Zf, x /Iq x pf,. 

Corollary 2.4. The lifted Block ind-variety B = IJ(fc^)ez2 ^k/ = Uk e e^k,e,Q 
is a reduced ajfine ind-variety and its components are nonsingular subvari- 
eties 

BkAo ■■= {(^1. ^2, A) G C* X C* X C I M^tn (P, ^) - A = 0} (15) 
where {m,n) E Za x md {pi, P2) G Pa x Pb- 

Obviously we have a group action of pa x Ph on each -Bfc/,(? by acting 
on each fiber p ■ (zi, ^2, A) = {piZi^ P2Z2, A), in other words the group action 
changes the powers q in Mmf^{p, z) where p= (pi, P2) = (Pod ^02)- 

2.3 Singularity locus 

In this subsection we describe explicitly the singularities of the components 
of the Bloch ind-varieties that we need to deal with. 

For fixed {k^tj and consider a typical fiber Ex = Ex{k,i,g) of -Ba;,^,^ 
given by the set 

{{z^, Z2)\X = p'^t^''^''''''zi+po"^^'''-''^''-^'z^'+pT^^^^^ 

(16) 

and take the derivatives formally 

dEx _ a{n+£b)a+p _ -a{n+eb)a-p -2 

~ POI POI ^1 

dEx _ P(m+ka)h+q _ -l3{m+ka)b-q -2 V V 

a57 ~ P02 P02 ^2 

where we use the same notation Ex for the variety and for the polynomial 

TP _ „0'in+^>')a+p -a{n+£b)a-p I B{m+ka)b+q -I3{m+ka)b-q I . 

-C/A - Poi + Poi ^1 + Po2 Z2 + Po2 Z2 - A 

Then the singular points consist of four points (ipo"*'"^^^^" ^, ±Po2'''''"'^'^"'''' '^). 

We have an analog of Lemma 5.1 of [5]. When A = 0, -Ea splits into two 
components 



il^l,^2j|Poi ^1+P02 ^2 - U| , s 

\\.Zl, Z2)\Pqi Zi+Pq2 Z2 - U| 
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When X = 4 , E\ is irreducible with (pg p-^i'^+^"-)^-i'^ 
only singular point. 

When \ = —A , Ex IS irreducible with (— _p-P'y'^+'^°-)^-i^ 
its singular point. 

Otherwise, when A G C \ {0, ±4}, the typical fiber is a nonsingular com- 
plex curve. 

2.4 Fermi curves 

In this subsection we identify the components of the Fermi curve associated 
to the Bloch ind-variety. It is these countably many components of the Fermi 
curve that will contribute to the period computation of §2.6[ after taking care 
of the compactification and singularities problem. 

Recall that affine Fermi curves are defined by -Fa(C) := 7r~^(A) with the 
projection vr : i? — )> C; (^1,^2, A) 1— )■ A. Since B = B / fj,a x fj,h, Fx{C) is given 

by 

[J {(^l,^2)|A = e^™('^+^'')^j^_)_e-2'r«a("+^'')^-l_|_g27rj/3(m+fca)^^_|_g-27ri/3(m+ 
k,£,m,n 

(19) 

where 1 < m < a,l < n < b and {k,i) running through 7?, so we have 
countably many components for each A. Then the Fermi curve itself 

is, in this case, an ind-variety. Namely, for fixed integers {k,i) and {m,n), 
the component p^'^'"^'^ is given by 

(20) 

and the singular locus of F^^' is easily derived from that of Ex- 
When A = 0, p^'^'"^'"^ has two components 

{(6,6)|e'"'"("+'')ei + e2-^''(™+'='^)6 = 0} . . 

{(6,6)|e^™("+^^^^i + e-2-*/^(™+'='^)^2"^ = 0} ^ ^ 

When A = 4, F^^'^'"^'" is singular only at (e-2™(«+^fe), e-2^*'^("^+*=")), simi- 
larly when A = -4, Fj^'^'"^'" is singular only at (_e-2™(«+^fe), _e-2^*/3(m+fca)) 
and they are irreducible curves. 
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2.5 Compactification and blowups 



In this more technical subsection we exphcitly describe the form of the com- 
pactification of the components of the Bloch ind- variety and the compatibihty 
of this compactification operation with the hmit procedure described in §2.11 
above. We also describe the blowup procedures that take care of the singular- 
ities problem and again check that these can be carried out compatibly with 
the operation of passing to the limit. Performing both of these operations 
will create a more complicated "double limit" procedure, which can be ap- 
propriately described, with all the compatibility conditions directly encoded, 
by the notion of an ind-pro system of varieties. The resulting double limit 
obtained in this way is called an ind-pro-variety and is then the geometric 
space that we need to deal with, which describes the complex energy-crystal 
momentum dispersion relation in the case of the Harper and almost Mathieu 
operators with irrational parameters. 

We proceed as in §4 of [5] to obtain the compactifications of the compo- 
nents of our Bloch ind-variety. We use the same notation and terminology 
as in §4 of [S]. 

By adding points at infinity, first we have Bi as the projective closure of 
S in X X Let s,t belong to the set {0, oo}, it is easy to see that 
the intersection Bi {^{F^ x P^ x P^ \ C* x C* x C) consists of eight rational 
curves {s} x P^ x {oo}, P^ x {t} x {oo}, {s} x {t} x P^ as in Lemma 4.1 of 
0. 

Bi is singular at the points Os,t '■= {s,t, oo), after blowing up these four 
points, we define B2 as the strict transform of Bi and the group action of 
fia X fJ^b can be lifted onto B2 naturally. 

Let P*'* be the exceptional divisor over Og^t, the intersection points of 
P*'* with the strict transforms of eight rational curves are denoted by Wq'* G 
P^'* n {s} X {t} X Pi \ 0,^t, wl'' G P^'* n {s} X Pi X {oo}\0,,t and w'/ G 
ps,t Pi pi X {t} X {00} \ Os^f Obviously *, wj'*, Wg'* are fixed points of 

We then have a direct analog of Lemma 4.2 of [5j adapted to our inductive 
system of varieties. The main difference lies in the fact that here we deal with 
a countable family of quadrics arising from the blowups instead of having just 
ab of them as in the original case of [5]. 

Lemma 2.5. There exist quadrics Qk\g C P*'* containing Wq^ , wj'*, 
satisfying that there are neighborhoods Ui of w^'^ such that Ui fl B2 consists 
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of countably many branches. 

Proof. We only consider what happens around the point (0, 0, oo), other cases 
can be treated similarly. 

Using the coordinate system Zi, Z2, ^ = in a neighborhood of (0, 0, 00), 
fix s = t = 0, then w^^ = (0, 0, 1), = (0, 1, 0), = (1, 0, 0). 

Recall that the closure of the lifted Bloch variety Bi is the zero-set of the 
determinant of the matrices fiziZ2{Mm^'n {p, z) — XI) with diagonal entries 

a{n+eb)a+p 2 ^ , -a{n+eb)a-p , I3(m+ka)b+q 2 -0[m+ka)b-q 

Poi A^^l ^2 + ^01 /^^2 + Po2 /^^1^2+P02 IJ,Zi — ZiZ2 

(22) 

Define 

'^k,l,Q ■— \\Zi . Z2 . H) ^ r I 11Z2 + Pq2 l^^l ~ ^1^2 — Uj- 

(23) 

therefore P™nS2 is the union of infinitely many quadrics Qk,e,Q- In addition, 
B2 is nonsingular at every point not being a point of intersection of different 
quadrics. 

Let us look at Wq^ first, from the tangent part "("+^'')'* ^^2+^02^*''"''''^"'''' 
there are infinitely many transversal branches labeled by {k,i), {m,n) and 
(p, q) intersecting at L :— {zi = Z2 = 0}. 

Then consider w^*^, since wf^ — (0, 1, 0) in the coordinates Zi,Z2, intro- 
duce u ^ 0,ui,U2 such that 

zi — uiu, Z2 — ly, n — U2iy 

And now B2 is given by the determinant of the infinite matrix with diagonal 
entries 

p-i-+a)a+p^2^^^2^p-ain+eb)a-p^^^^^^^^ 

(24) 

So we get a singular line L' := {ui = U2 = 0} and the tangent cone of B2 
at L' is the union of infinitely many planes {Pqi^^^^^^"'~^U2 — Ui = 0}. 
Blow up the line L' by defining new coordinates 

Ui = V1V2, U2 = V2 

The strict transform of B2 then is given by the determinant of the infinite 
matrix with diagonal entries 

a(n+£b)a+p 2 2 2, —a(n+(.b)a—p , 8(m+ka)b+q 2 , —3(m+ka)b—q 

POI V^V2U + Poi +P02 V1V2U + V1V2 - Vi 

(25) 
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In other words, 



-a{n+£b)a-p -I3(m+ka)b-q -a{n+eb)a~p 

POI - ^^1 + Po2 Poi ^2 /r,px 

-B(m+ka)b-q / -a(n+eb)a-p \ , i • i V^"/ 

— Pq2 iPol ' — Vi)V2 + higher corrections 

We get countably many singular hnes Li^n,p '■= {vi = p^^^"'^^^^^ ^ ■,V2 = 0} 

1,1, , J. 1 r • / —a(n+£b)a—p \ , -~B{m+ka)b—q —a{n+lb)a- 

and the tangent cone at each L£ ,^ p IS (Pq]^ '^ij+Po2 Poi 

i.e. it has countably many branches as claimed. 

Finally consider = (1,0,0), analyze similarly to wf^ and the strict 
transform B2 is determined by entries 



-(3{m+ka)b-q -a{n+lb)a-p -I3(rn+ka)b-q 

P02 —Vl+ Pg2 f 2 

—aln+lb^a—p / —B(m+ka)b—q \ , i • i 

— Poi (,Po2 ~ ^1^2 + higher corrections 



(27) 



Countably many singular lines Lkrn,q '■= {vi = Pq2^^^''"'^^ '^,V2 = 0} are 
obtained in this well. Let denote these singular hnes Li^n,p and 

Lk,m,q together. 

□ 

By Bezout's theorem Q^' ^ ^ and Qlfii^^i have a point of intersection rf*'* 
other than Wq'*, wJ'*, Wg'* if and only if {£,n,p) 7^ {£',n',p') and {k,m,q) 7^ 
(fc', m', q'). If such c?'^'* does exist, then it can be proved that they are ordinary 
double points of B2 as in [5]. 

As in the above lemma (12.51) . we blow up the strict transforms of the 
rational curves {s} x {t} x P\ {s} x x {00} and x {t} x {00}, and let B3 
be the strict transform of B2. Then B^ is nonsingular on the strict transforms 
of the rational curves {s} x {t} x P^, since all branches are transversal to 
each other in this case. On the other hand, it is singular at countably many 
lines which lie over the strict transforms of {e} x P^ x {00} and respectively 
of P^ X {/} X {00}. We now take care of these singular lines Lk/^g. 

In order to do this, we first need to clarify how we are going to proceed at 
performing the "double limit" that will result when we keep adding more and 
more components to the Bloch variety (constructing the ind-variety) and at 
the same time blowing up all these singular lines, in a compatible way. The 
double limit procedure we need is encoded in the notion of an ind-pro-system 
of algebraic varieties. Again, we first recall the general abstract definition of 
this procedure, then we give a more heuristic interpretation of its meaning, 
and then we apply it to our concrete case to obtain explicitly the construction 
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of the compactification and resolution of the Bloch variety in the irrational 
case. 

Definition 2.6. An ind-pro system of varieties {Xf } is a double indexed 
set of varieties such that for each i > 0, {Xf } is an inductive system while 
for each n > 0, it is a projective system and every square of the system is 
Cartesian, i.e. the horizontal maps are injections and the vertical ones are 
surjections, the diagram commutes and that the top-left corner is the fibered 
product of the bottom and right map. 

Xmc 

Xf^ >XJ 

Using categorical dual notions, we define a pro-object of the category 
Ind(Var) as a formal cofiltered limit of a projective system of ind-varieties. 
Further, we can construct the category ProInd( V^ar), or ProKoIndKo(^ar) if 
the index sets are countable. 

Given an ind-pro system of varieties {X^}, taking the inductive limit and 
the projective limit gives the ind-pro- variety ^m . li^^ X", which is an object 
in ProKoIndKo (V^ar). 

What this abstract definition is saying in more heuristic terms is that 
we are building up a space by performing simultaneously two kinds of op- 
erations on a collection of algebraic varieties. One type of operation is the 
one we have already seen when discussing the ind- limit, namely we orga- 
nize our varieties in a sequence of inclusions that progressively adds more 
and more components. The other operation takes care of the blowups that 
progressively remove singularities (these form a projective system, since the 
smoother blowups map down by projections to the more singular varieties 
they are obtained from) . The problem here lies in the fact that the required 
sequence of blowups needs to be performed on an infinite number of com- 
ponents that intersect each other, and the way to make this possible is by 
compatibly carrying out blowups on the approximating varieties of the induc- 
tive system. Notice that, to make sense of this double limit, we only need 
to use the very simple categorical notions of inductive and projective limit. 
In fact, the categorical definition above is simply encoding the compatibility 
condition that makes this kind of double limit operation possible: it is say- 
ing that operations we perform at the level of the approximating algebraic 
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varieties X", as long as they are done in a way that is compatible with the 
maps between the varieties that define the two limit operations, will carry 
over to the limit, even though the space we obtain as the double limit of this 
sequence of algebraic varieties is no longer an algebraic variety but a Cantor- 
like geometry, which cannot be directly described within classical algebraic 
geometry. 

We now return to our specific case of the Bloch varieties for the Harper 
operator in the irrational case. 

Proposition 2.7. In the case of irrational parameters, the compactification 
of the lifted Bloch ind-variety, is an ind-pro- object defined by a chain of 
iterated blowups on the ind-variety B. 

Proof. For fixed {k,i),{m,n) and {p,q), we can blow up the line or 
Lk.m,q as usual. Let us look at I/^,„,p here. By introducing new coordinates 
Wi, W2 such that 

—a(n+£b)a—p /r>o\ 
^^1 - POI - ^1^2, V2 = W2 (28) 

the defining equation becomes 

-a{n+eb)a-p -I3(m+ka)b-q , -a{n+B)a-p I3{m+ka)b+q 2 , -0{m+ka)b-q 

Poi P02 — Wi-\- Poi P02 + P02 W1W2 

^^-ain+eb)a-p^^^2 ^ 2wiwiiy'' + p^^^^''^'^'' W^W2V^ + p"^'^^'''^''^'' ^^2^^ 

(29) 

First fix {k,(!^) = (0,0) and let g = {m,n,p,q) vary, we get 2ab singular 
lines {LQ^n,p, Lo,m,q} , then blow up these 2ab lines, define as the strict 
transform of B^. We start over this process again, now let \k\ < 1, \i\ < 1 and 
g vary, blow up these Gab singular lines IJ|fc|<i \(:\<i{^k,i,s}, define Ai as the 
strict transform of B-^. In general, we have Ai being the strict transform of B^ 
after blowing up 2(2i-|- l)a6 singular lines U|jk|<i \e\<i{-^k,i,e} for |A;| < i, \i\ < i 
and all g. 

Obviously, the projection maps tt^ : Ai — > Aj are surjective for any 

^ ^ J ^ 0. In other words, we have a projective system of ind- varieties since 
each Ai is an ind-variety derived from the blowup a : Ai ^ B^ and -B3 being 
an ind-variety. Hence, the projective limit A = ]miAi is an ind-pro- variety 
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in ProKoIndKo(V^ct'")- We have the following commutative diagram 



A 




and it is easy to see our construction satisfies the Cartesian squares. □ 

We call A the compactification of the lifted Bloch ind-variety B and 
denote it by B from now on. 

What this result shows is that we can compatibly resolve the compactifi- 
cation and singularities problem on each individual finite approximation to 
the Bloch ind-variety, in a way that allows us to pass to the limit. Now, 
as a consequence of this construction, we find that the space we obtained 
as the double limit of this compatible family of algebraic varieties is indeed 
a Cantor-like space, but one that has a good approximation by algebraic 
varieties. In fact, Kapranov showed in [9] (see also [15|) that the category 
of locally compact Hausdorff totally disconnected spaces can be identified 
with a full subcategory of IndPro(S'eto), which means that geometric spaces 
obtained by the kind of double limit procedure described above look Cantor- 
like. One can find a sketch of the proof in [15]. In particular, this means 
that the compactification of the Bloch variety of the Harper operator with 
irrational parameters is an ind-pro- variety, which is a Cantor-like geometric 
space, as one might have expected by thinking of it as the complex energy- 
crystal momentum dispersion relation, in a case where the band structure in 
the spectrum is replaced by a Cantor set. The important additional informa- 
tion that the result above gives us is the fact that this Cantor-like geometry 
admits a good approximation by algebraic varieties: this will be useful in 
§2.61 and the rest of the paper below, since it will allow us to reduce the cal- 
culation of the density of states for this Cantor-like geometry to a sequence 
of terms that can be computed compatibly over the approximating algebraic 
varieties, and that can therefore be identified explicitly with period integrals. 

We now rephrase the previous result on the Bloch varieties in terms of 
Fermi curves, since these will be the curves over which the period calculation 
of §2.61 for the density of states will take place. 
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We have the analog of Theorem 4.2 of [3]. The map vr o c : B — t- C 
extends to a morphism vr : i? — )■ and its fibers Fx '■= 7f^^(A) are called the 
lifted Fermi curves. The inclusion map i : B "-^ C* x C* x C gives rise to a 
morphism z : S — )■ x x 

Theorem 2.8. B \ B consists of countable curves 
(i) sections ofii, S^'^^ with ?(S^'^^) = {s} x {t} x P^ 

(a) quadrics Ql'),^^ with iiQ^^f^^) = {s} x {t} x {oo} 

Hln,p with i{Hl^ J = {s} X pi X {oo} and if* „,,^ with i{Hl^J = 
P^ X {tj X {oo} 

These curves meet transver sally at only one point, except that {rf*'*} are 
ordinary double points of B, every intersection point is a nonsingular point 
of the compactification B . 

Taking the quotient by the structure group /Iq x /i^, we get the compact- 
ification B = B/ fia X fJ'b of the Bloch ind- variety B. At the same time, we 
get two morphisms, the inclusion z : 5 — )■ P^ x P^ x P^ and the projection 
TT : 5 — )■ P^, whose fibers are the so-call compactified Fermi curves Fx- Let 
^M,m,ri' Qk/,m,n^ ^^^1 Hj,^^ in B be the image under the quotient map, 
in addition, let c?*'* be the image of c?'^'* in B. 

Theorem 2.9. B\B is the union of countable curves, 

(i) ^kAm,n '^^^^ '^i^ki,m,n) = {^} ^ {0 ^ ^'^ > nonsingular points of B 

(a) Q'kimn with ordinary double points at (i'^'* and nonsingular at all points 

of Qt^Ay^d''} 

(ill) Hl^ with i{HlJ = {s} X Pi X {00} and with i{HlJ = ¥'x {t} x 
{00}, nonsingular points of B. 

2.6 Density of states and periods 

In the above subsections we have constructed the desired geometric space, 
the compactification of the Bloch ind-variety, which can be approximated by 
finite dimensional algebraic varieties. In this subsection, we will proceed to 
compute the density of states on those approximating components for the 
Harper operator. The density of states is indeed the period integral over 
Fermi curves, which also related to the periods of elliptic curves. 
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In this subsection, we assume |^i| = |^2| = 1, i-e. (65^2) = (e^'^*'^!, e^''*'^^) 
for some {ki, ^2) G (0, 1]^, then the Bloch ind-variety is given by 

%l){m + a,n) = e^'"'^^^ip{m, n), (30) 
?/'(m, n + h) = e'^'^^''^ip{m, n)} 

Denote its spectrum by cr{H) := {Ej{ki, ^2), j G N}, the function Ej{ki, /C2) 
is the so-called j-th band function. 

In order to compute the density of states, we get back to a continuous 
model by taking the limit of the lattice model. Here again we use the same 
notation and terminology as §3 and §11 of [5j. 

Let Hn denote the Harper operator H acting on ['^{7? /ariL © 6nZ) for 
some integer n > 1, later we will take the limit as n tends to infinity. It is 
easy to see that the eigenvalues of Hn are just given by 

mi m2 .. 

{Ej[ — , — ) 1 < mi,m2 <n, j >l \ 
n n 

To define the integrated density of states, we first count the number of 
eigenvalues less than or equal to A. Or equivalently, with the step function 
9(x), define 

00 n 

-n{X):=Y: E 0(A-i?.(^,^)) (31) 

j=l mi,m2=l 

The integrated density of states is then defined as the limit 

p(A) := lim -^i^niX) (32) 

and the density of states is defined as the derivative dp/dX. In the literature, 
the integrated density of states can also be defined as the normalized trace 
on a III factor in von Neumann algebra theory [18]. 

Rewrite it as an integration, we get an integral over the continuous vari- 
able k e [0, i]2 



p(A) = lim„^^ J, Er=i ^ T.l,,m.=i 0(A - E,{^, ^)) 



(33) 



Lemma 2.10. The density of states can be expressed o-s ^ = ^ Ix&a{H)^^' 
where uj\ is a differential 1-form. 
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Proof. 



1 voo r 



ab j=lJEj{k)=X IVfcEjl 



J^VOO r y/dkf+dk^ 



_ yoo r 

ab j=l JEj(k)=X ^diE]+d2E'f 

For fixed A, from Ej(ki, /C2) = A, we have dEj = ^-dki + ^^dk2 = 0, 
dki = — ^ — -dko or dko — — ^ J dki 

Then 

ds = ^d,E^ + d^E] ^ or ds = ^d,E^ + d^E] ^ 
Therefore we can write the density of states as 

dp 1 /" dk2 1 /" dk\ 1 /" (ifci 



d>^ ab jr^i J E^{k)=x \9iEj\ ab J E.{k)=x \92Ej\ JAeu^{E^(fe)} 

Denote the differential 1-form by ux, i.e. u}\{ki,k2) := = \d2E \ ' 

lemma is proved. □ 

Let P(^i,,^2, A) be a general polynomial with ^1 = e^'^*'^^ and ^2 = e^'^*'^^. 
Then d/ci = d^i/2-Ki^i and dfc2 = d^2/'27ii^2, we also have dA = diEjdki + 
d2Ejdk2 from A = /cs). Plug into dP = P^.d^i + P^^d^2 + PxdX = 0, 

we get 

(2mP^,^i + PxdiEj) dki + (27riPj,6 + Pa52^j) dk2 = (34) 
Since dki and (iA;2 are independent, so 

27riPai + PxchE,=0 

2mP^,^2 + Pxd2Ej^0 ^''''^ 
Lemma 2.11. The density of states is a period integral over Fermi curves. 

Proof. Recall that the Bloch ind- variety B = U(^k/)&z'2 ^'^n=i Bm,n , where 

^{k,t} ^g27ria(n+^&)^^_j_g— 27ria(n+€6)^— l_j_g27ri/3(m+fca)^^_j_g— 27ri/3(m+A;a)^— 1 ^ gj^ 

(36) 
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and Fermi curves Fx = U(fc_^)ga2 U'^^^^i pk,t,m,n^ where 

rpk/,m,n f 2-Kia{n+lb) (: , —2iTia{n+eb) 2iTil3{m+ka) , —2iTil3{m+ka) \i 

-^A ~ 1^ Sl"rC -t-e C,2-r(i ^2 ~ ^/ 

(37) 

pk,£,m,n ^2nia{n+£b) _j_ ^—2nia{n+£b)^—l _j_ ^2nil3{m+ka) _|_ g— 27ri/3(m+fca)^— 1 ^ 

denote the polynomial, then 

pk,£,m,n 

pk,e,m,n _ ^2'rria{n+eb) _ g-27ria(n+^fe)^-2 {'^8) 
pk,e,m,n _ 2ni(^{m+ka) _ -27ri/3(m+A;a) t-2 

It follows that on each component i^^'^'™-'"- 



dT-i o • r>k,£,m,n ^ 
lEj = 2mP,l ^1 

Dk,{ 

6 



2Ej = 2mP,' ^2 



(39) 



Furthermore 



dk2 _ d^2/2TTi^2 _ _J- dg2 

diEj 27riP^'='''™'"a (27ri)2 ^^^^pk,e,^,n 

dk-t ^ <igl/27rtgi ^ dji 1'*^) 

i.e. wa(6,6) = ^^2|^^g^pfc,f,,n,n| = ^^2|g^^^p\A»n,n| - Bj the change of variables 

(^1 = e^''*''^ ^2 = e^''*''^, the graph (fci, A;2, Ej{ki, ^2)) is changed into (,^1, ,^2, A). 
Then 

^ = 4-^/ WA(A;i,fc2) = ^ y2 yi [ ^a(^i,6) = A;/ ^^a 
dA ah ^ J p.(u\_x ah ^ ^ ^ J pk,i,m,n ah Jp^ 

□ 

Notice how in this result we used explicitly the fact that the geomet- 
ric space describing the complex energy- crystal momentum relation for the 
Harper operator with irrational parameters admits a good approximation 
by a family of algebraic varieties, as we proved in §2.H - §2.5[ In particular, 
we see here that the density of states is computed as a period integral with 
compatible contributions from each of the components i?^'^'*"'" of the Fermi 
curve, as described in Theorem 12.81 above. 
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Recall then that the incomplete elliptic integral of the first kind is defined 

as 

F{x,k):^ = (41) 

and the complete elliptic integral of the first kind is defined as 

K{k) .-F{l,k)- -j======= (42) 

The complete elliptic integral of the first kind is also called the quarter period. 
We rewrite ujx on F^^'^''"'" n {1^1 = |6I = 1} as 

d£i dki 



47r2|^i^2(e2'^^'^("*+*^«) - e-2'^^/^("^+H^j2^ | 47r|sm27r(A;2 + ^(m + ka)) \ 

(43) 

Let C^'n^ {{ki,k2)\cos2-n{ki + a{n + ^h)) + cos2T:{k2 + I3{m + ka)) = A/2} 
so that 

^ = J_ V V / (44) 

47ra6^;^^^^^yci*/) |sm27r(A;2 + /3(m + H 

Theorem 2.12. If we define the elliptic modulus as k := |^|^; then the den- 
sity of states as a function ofk is ^\p,kAm,n = j^{l+k)K{k) = 2^-^(f+i) 
on each component of the Fermi curve. 

Proof. Denote the definite integral I^i\X) := /^(.,,) |,,„2.(fc,S(m+fca))| 

/ci e / moving on the curve Cm',n ■ 

The differential form cja is well defined, and ki ^ —a{n + £b), —a{n + 
ib) ±1/2 {mod 2n), we can just assume < —a{n + ib) < 1/2. 

Let ki run through a half period of the sine function, i.e. 27r(A;i + a{n + 
ib)) e (0,7r), then Jm,n''(A) is twice of this integral over a half period. 



^ ' ^l-(A/2-cos27r(fei+a(n+^6)))2 

_ Q rl/2 dcos27r(A:i+a(n+^&)) 

- ^ Jo 



Y^l-(A/2-cos27r(A;i+a(n+^b)))2(-27rsm27r(fci+a{n+^fe))) 
-1 rl/2 dcos27r(fci+a(n+^fe)) 

~ Jo 



^Jl-{\/2-cos2■n{k^+a(n+^h)))'^^Jl-(cos2■K{k^+a{n+a))Y 

dx\ 

Vl-(A/2-a;i)2^/l^ 
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in the last line we set Xi :— cos2Tr{ki + a{n + £b)), here —l<Xi<l and 

— 1 < — A/2 < 1, there are two cases: 
(i) A/2 - 1 < xi < 1 when A G (0,4); 
(u) - 1 < xi < 1 + A/2 when A G (-4, 0). 
In the first case, A e (0,4) 



7-(fc/)/\\ _ 1 fl dxi 



X 



V4-1 ^l-(A/4-x)2^1-(a;+A/4)2 
1 rl-A/4 

7^ Ja/4-1 ^((i_a/4)2-s2)((1+A/4)2-x2) 
1_ rl-A/4 dx 

■K J\ 



'A/4-1 (i_a/4)(1+A/4)^(1-x2/(1-A/4)2)(1-x2/(1+A/4)2) 

i:=x/(l-A/4) = Li ^(i_,2)Ji_fc2,2) 

A;:=(4-A)/(4 + A) = i±^i^(A;) 

In the second case, A G (—4, 0) 



^l-(A/2-a:i)20^ 

x:=x,-A/4 



t:=x/(l + A/4) =n^/_\ 



^l-(A/4-a:)2^1-(a;+A/4)2 



(4_A)7r J-1 ^(I_t2)(i_jfc,2i2) 

A;' := (4 + A) /(4 - A) = ^K{k') 

By identifying k — k', namely, defining k := we conclude that 

/^;^)(A) = ^i^(A;) and the density of states on each component F^'^'""'"" 
turns out to be 

^K(k) - -^K(^) 

□ 

Interchanging positive and negative A corresponds to the electron-hole 
symmetry in solid state physics. We give a representation of the density of 
states in terms of elliptic curves. 

Corollary 2.13. On each component of the Fermi curve the density of states 
IS a sum of two half-periods of isomorphic elliptic curves up to a constant, 
namely = ^^U^^^k + /y ^i/fc)- 
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Proof. First we rewrite the density of states on each component as 



47r^a6 J-fe ^(jfc2_„2)(i_„2) 47r^a6 J-1 ^(I_t2)(i_fe2t2)/jt2 



ATr'^ah J—k , /(i-2__-\\(f2 _i.2\ in'^ah J — 1 



dt 



For the eUiptic modulus A; e (0, 1), consider the eUiptic curves given by 

Ek ■.= {{t,y)\y' = {t'-l){t'-k')} 
Eyk ■.^{it,y)\y'^ie-l)ie-l/k')} 

If we denote the canonical holomorphic form u = dt/y hy cuk for (resp. 
oji/k for Ei/k), then 



dp. 1 



Recall that for E^, we cut the Riemann sphere from A; to 1 and from 
— 1 to — /c, then assemble a Riemann surface homeomorphic to the torus T^. 
Similarly for -Ei/fe, we cut the Riemann sphere from 1 to 1/k and from —1/k 
to —1, also assemble a Riemann surface homeomorphic to the torus T^. 

A framed elliptic curve (i?, 5, 7) is an elliptic curve with an integral basis 
for the first homology such that the intersection number 5 ■ ^ — 1. And for 
the period vector ( cj, a;), the ratio t{E, 5, 7) = j^uj/ uj is an invariant 
for isomorphic complex tori, which is called the modulus of the isomorphism 
class. 

On E^, we can choose 5 as the cycle from A; to 1 then from 1 back to k 
and 7 from —k to k then back to —k, which makes {E^, S, 7) a framed elliptic 
cTuvc. Similarly for Ei/k, we choose S' as the cycle from 1 to 1/k then back 
to 1 and 7' from —1 to 1 then back to —1. 

By deforming the path of integration, Uk = J^^ cok + ^^k = ^k 
and J^ujk = j'lk'^k + jk^ ^k = '^ j-k^k- Then 
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It is easy to see that T{Ek) = r(_E'i/fc), indeed J_^ui/k = k J_i^ujk and 

Ji ^i/k = ^ Xc ^k- In other words, we have elliptic curves Ek isomorphic 
to -Ei/fc and we can interchange between them by rescaling the fundamental 
domains. 

Hence on each component of the Fermi curve 

□ 



2.7 Spectral functions 

With the density of states in hand, we can continue to derive some interesting 
spectral functions. In this subsection, we will calculate the partition function 
of the propagating electron in magnetic field on the components of the Bloch 
ind-variety. 

Fix {k,i,m,n), we consider one irreducible component Bm',n , other com- 
ponents can be treated similarly. Recall that 

TD{k/) f 2-Kia(n+lb) /: . —2nia{n+£b) . 2T7il3{m+ka) c . —2ni(3{m+ka) \ r\\ 

^m,n ~ Sl"rC "he C,2~r(i C,2 A — U| 

(45) 

For the moment, we omit the superscripts and denote the above defining 
polynomial as P(^i,^2,A). 

From dP = P^^ci^i + P^2dC,2 + P\d\ = 0, then on i?m,'n the pullback 

*rj\\ P^i^^^l P^2(^^2 r> It , r> Jt 
^ [dV = -5 -5 — = P^id^i + ^6^6 

If we wedge this form with 

1 



i2m)Hi^2P^2 
then 

n.A.-(dA)^ J '^'f (46) 

This observation was already made in §11 of [5j. Note that fix is slightly 
different from ux. As defined before, ux is a volume form so it should be 
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positive and we finally got a positive period. By contrast, as for Q,x, we get 
rid of the absolute value and take the orientation into account. 

Define Q := QxA7r*dX over Bm',n such that Q\pk/,m,n = Q\, where /t"^'^'™'" 

is the component of the Fermi curve of Bm'n, ^ is called the density of states 
form and then Qx is the relative differential form with respect to 7r*dX. 

^ /o<|A|<4 J^fcA-." '^A 
~ 1^ Jo 1+k "'^ 

Then it is possible to construct some interesting spectral functions based 
on dp for each component. Let us look at the zeta function of the Harper 
operator on Bmln, 



[ X'dp= I (47) 

Jo<|A|<4 Jb^'^J 



/o<|A|<4 ^^'^P - lo ^^d^ Jpk,e,m,n UX + X^dX JpkAm,n Ux 

~ TT^ab Jo\l+k) 7r2a6 JoVl+fc^ 1+fc "'^ 

In particular, when s = 2/c + 1 is an odd integer, ('^^'"^'"'(2k + 1) = and 
when . = 2fc is even, Cj'--(2fc) = |^ I^iUf'^dk. 

There is another way to compute the zeta function using two dimensional 
residue theorem. Denote 

namely, P(^i, ^2, A) = VF(^i, ^2) — A. Since for any function /(^i, ^2, A), 
/t^o<|a|<4 6, A)5(A - W{^,, ^2mid^2dX 

^ /t2,0<|A|<4,A=w/(^1'^2, A)|^ 
~ Jt2,0<|A|<4,A=W 'S2> /l+p2 ^p2 
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Therefore 



By the residue theorem, only the positive integer powers survive and the con- 
tour integral only depends on terms with 1, ^1, ^2, and ^1^2 in the integrand. 



We can also have the partition function of the Harper operator on Bm, 



n , 

r -.00 ^k,£.m,n / r\-r,\ 

In order to get the spectral functions on the whole variety we should 
collect all the contributions from countable components. 



3 Almost Mathieu Operator 

We compare in this section the case of the density of states of the two- 
dimensional Harper operator analyzed above with the analogous problem for 
the one-dimensional almost Mathieu operator. In particular, we show how 
to recover in the one- dimensional case the familiar picture of the Hofstadter 
butterfly and the corresponding density of states with its explicit dependence 
on the parameter. 

3.1 Algebro-geometric model 

We apply the same process to the almost Mathieu operator, the big difference 
is that now its density of states depends on the parameter a. Therefore the 
derived spectral functions will have totally different properties compared with 
those of the Harper operator. Here again we have the analog of §3 of 0. 
First the Bloch variety of the almost Mathieu operator is now given by 

B' := {(e, A) G C* X C I H'^in) = Xip{n), ip{n + a) = ^ip{n)} (49) 

It is easy to see that (^, A) belongs to B' if and only if A is an eigenvalue of 
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the matrix M' with countable components M'^, £ G Z 








V 



-1 







1 




In other words, 



= {(e, A) G C* X C I det{M' - A/) = JJ det{M'^ - A J) = 0} (51) 



In fact we can expand the determinant as det{M'^ — AJ) = p£(A, a) —^ — ^ 



Pi{X,a) := + {2j2cos27ra{j + ia)){-\y'^ + OiX"-^) (52) 



Consider a continued fraction expansion {c„ = — } to approximate the 
irrational a. Then we let the prime a vary according to {qn} so that the 
spectrum of the Bloch variety can be approximated by numerical computa- 
tion of finite band structures. Indeed, based on semi-classical analysis and 
renormalization method, the Hofstadter butterfly is obtained by numerical 
approximations p]. 

Define the unramified covering and B' := c'^^{B') as before 



The structure group fia of the covering d : B' ^ B' acts on the fibers as 



Similarly, we use Fourier transform to convert the spectrum problem into the 
study of the Fourier modes, then the almost Mathieu operator is represented 
by the infinite matrix 

M' = Diagipz + p-^z'^ + 2cos2-Ka{] + la)), 1 < j < a, £ G Z (53) 

So under this representation, the lifted Bloch variety 



B' = {{z, A) I n n n^-^^ + P^'^^' + 2cos(27ra(j + £a)) - A) = 0} (54) 



a 



c' : C* X C ^ C* X C 
(z,A) ^ (^»,A) 



P- {z,\) = ipz,\) 



a 



i& pei^a j=l 
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and the Bloch variety 

a 

B' = {(e, A) I n + + 2cos27ra(j + £a) - A) = 0} (55) 

^ez i=i 

In this case, the Fermi curves are degenerate points 

a 

= U I ^ + + 2cos27ra(j + ia) = A} (56) 

3.2 Density of States 

Let us look at the self-adjoint boundary value problem in dimension one. For 
k & I, we now consider 

B' := {(e2'^'^ A) | H'ip{n) = \ip{n), ip{n + a) = e^^'^^{n)} (57) 

Assume its band functions consist of {Ei{k)}. 

Let H'^{n > 1) denote the self- adjoint operator H' acting on /^(Z/anZ). 
For k G ^Z, Ei{k) is an eigenvalue of H'^ and the spectrum of H'^ is 

TTX 

{Eii—)\i > 1, 1 < m < n}. 
n 



Similarly we have 

oo n 

,m 
n 



^n(A)-EEQ(^-^^(-)) 

Then the integrated density of states of the almost Mathieu operator is 
p'{X) = lim„^^ i ET=i El=, 0(A - E,{^)) 

= \YZi!M>^-m))dk 

And the density of states of the almost Mathieu operator is 



(58) 



l-E/^(^-^.w)*-E^ (59) 

where gi{k) := A — E^ik) and ki are the real zeros of gi{k). From the Fermi 
curve of the almost Mathieu operator, we have Ei{k) = Ej{k) for some j, i 

E\ (k) = 2cos2TTk + 2cos2TTa{j + £a) = A (60) 
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i.e. ki G I satisfies the equation 2cos2Tcki + 2cos2TTa{j + ia) = A. It is easy 
to see that 

g'.{k) = -—E^.{k) = A'Ksin2'Kk 

CLK 

Hence 

dp' 1 1 1 V ^ 1 

d>^~ 47r|sm27r/ci| ~ 47ra ^ ^ y^l - (A/2 - cos2na{j + ia)f 

Therefore the density of states of the almost Mathieu operator is a func- 
tion of the parameter a. However, as we have seen that the density of states 
of the Harper operator is independent of the parameters a and /3. 

In dimension two, the integral variable was absorbed by integrating over 
one- dimensional Fermi curves and the dependence of the parameters was 
resolved by the symmetric form of the magnetic translations. But for the 
almost Mathieu operator, the Fermi curves are degenerate points and the 
associated measure is just the counting measure, so the dependence of the 
parameter still remains, as expected, in agreement with the form of the den- 
sity described for instance in [7]. 



3.3 Spectral functions 

Now we can similarly define the zeta function of the almost Mathieu operator 
on each component B'^j 

(61) 

By the residue theorem, only the positive integer powers survive and the 
contour integral only depends on terms with 1 and ^. So consider powers 
n G N and expand the polynomial as 

{i + + 2cos2'na{] + £a))" = ^ C'^'^^^^^^i^^-^^{2cos2-Ka{] + la)f' 
There are only two cases 

(l)A;i — A;2 = and ki + k2 + k3 = n {2)ki — k2 = l and ki + k2 + k3 = n 
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thus Cwin) can be computed as 'winding numbers' 

- 2^ (k\Y(n-2k)\ ^ klik + l)\{n - 2k - l)\ 

n-2fc>0 \ ' \ ' n-2fc-l>0 \ J \ I 

To get the zeta function of the whole variety B\ we have C^H'ip) = 
X]^ez X]j=i Ch'' ('^)- Finally, we also have the formal partition function of 
the almost Mathieu operator on B' 



k=0 



4 Conclusions 



In the case of a discretized periodic Schrodinger operator describing electron 
propagation in solids, the complex energy-crystal momentum dispersion re- 
lation is described geometrically by an algebraic variety, the Bloch variety, 
which consists of the set of complex points that can be reached by analytic 
continuation of the band functions. The density of states can then be com- 
puted p] as a period on a curve (the Fermi curve) in the Bloch variety. 

In this paper we consider the case with magnetic field, where the periodic 
Schrodinger operator is replaced by the two-dimensional Harper operator, or 
its degeneration, the one-dimensional almost Mathieu operator. One knows 
from the spectral theory of these operators ([7], [13], [18]) that, in the case 
of irrational parameters, the band structure of the spectrum is replaced by a 
Cantor set, giving rise to the well known Hofstadter butterfly picture. Thus, 
one can see that, correspondingly, the geometric locus describing the com- 
plex energy-crystal momentum dispersion relation and replacing the Bloch 
variety will no longer be directly described by classical algebraic geometry. 
However, as one approximates the Cantor-like spectrum by a family of inter- 
vals (for example by approximating the irrational parameter by rationals via 
the continued fraction algorithm), it should be possible to correspondingly 
"approximate" this geometric space by ordinary algebraic varieties. The first 
part of this paper consists of a geometric result, which shows exactly what 
this approximation and limit procedure consists of. In particular, we show 
that one can obtain the space describing the complex energy-crystal momen- 
tum dispersion relation for the Harper operator with irrational parameters 
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as a "double limit" of a family of algebraic varieties, or more precisely as an 
ind-pro-variety, where one limit takes care of the presence of infinitely many 
components and the other limit of the blowups that are needed to deal with 
the singularities. The resulting space has indeed a Cantor-hke geometry, but 
one that admits a good approximation by algebraic varieties, in the sense of 
this double limit procedure, so that one can still use methods from classical 
algebraic geometry, applied compatibly to the varieties in the approximating 
family. 

We then use this geometric result to show that we can still obtain an 
explicit calculation of the density of states for the two-dimensional Harper 
operator as a period on the Fermi curve, where the period integral now 
consists of a sequence of compatible contributions from the components of the 
approximating system of varieties. This integral is then explicitly computed 
in terms of elliptic integrals and periods of elliptic curves. Similarly, we obtain 
explicit formulae for the spectral functions, again in terms of compatible 
contributions from the approximating family of algebraic varieties. We apply 
the same technique to the density of states and spectral functions in the 
case of the one-dimensional almost Mathieu operator. The main difference 
between the two cases is that, in the two-dimensional case the dependence on 
the parameter disappears in the density of states, because it is absorbed in the 
integration over the Fermi curve, while in the one-dimensional case the Fermi 
curves are points with the counting measure and one recovers the density of 
states obtained, by different methods, in [7], with its explicit dependence on 
the parameter. By a residue calculation, we also obtain the zeta function of 
the almost Mathieu operator as a sum over "winding numbers" associated 
to the components in the approximating family of algebraic Bloch varieties. 
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